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ABSTRACT 

We propose a strategy to measure the dark matter power spectrum using 
minimal assumptions about the galaxy distribution and the galaxy-dark matter 
cross-correlations. We argue that on large scales the central limit theorem 
generically assures Gaussianity of each smoothed density field, but not coherence. 
Asymptotically, the only surviving parameters on a given scale are galaxy 
variance a, bias b = Q s //3 and the galaxy-dark matter correlation coefficient r. 
These can all be determined by measuring the quadrupole and octupole velocity 
distortions in the power spectrum. Measuring them simultaneously may restore 
consistency between all /3 determinations independent of galaxy type. 

The leading deviations from Gaussianity are conveniently parameterized 
by an Edgeworth expansion. In the mildly non-linear regime, two additional 
parameters describe the full picture: the skewness parameter s and non-linear 
bias 62 • They can both be determined from the measured skewness combined 
with second order perturbation theory or from an N-body simulation. By 
measuring the redshift distortion of the skewness, one can measure the density 
parameter Q with minimal assumptions about the galaxy formation process. 
This formalism also provides a convenient parametrization to quantify statistical 
galaxy formation properties. 



1. Introduction 

The measurement of the distribution of matter in the universe has been one of the 
frontier goals of modern cosmology. The correlation of galaxies has been measured in 
several surveys and is known to significant accuracy. The abundance of data has led to 
many theoretical challenges, especially for flat Cold Dark Matter (CDM) cosmologies. In 
the simplest models, galaxies are considered an inbiased tracer of the mass. Several different 
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measurements of velocities then allow us to measure the density of the total matter, 
using either cluster mass-to-light ratios or pairwise velocities (Peebles 1993). Under these 
assumptions, one obtains values of the density parameter Q ~ 0.2. It is known, however, 
that galaxies of different types correlate in different ways (Strauss and Willick 1995). 
Thus, not all galaxies can simultaneously trace the mass, and it appears plausible that no 
galaxy type is a perfect tracer of mass. At this point, galaxy formation is not completely 
understood. Simulations of galaxy formation suggest that the galaxy formation process is 
stochastic and non-linear (Cen and Ostriker 1992). This complicates the derivation of f2 
from dynamical measurements. Peculiar velocity fields in principle allow one to measure 
the mass fluctuation spectrum (Strauss and Willick 1995), but even here one must assume 
that statistical properties of galaxies, such as the Tully-Fischer relation, do not depend on 
the local density of matter. 

Assumptions need to be applied to relate the distribution of visible galaxies to that of 
the total matter. The most popular of relations has been to assume linear biasing, where 
the power spectrum of galaxies P g (k) is related to the matter power spectrum P(k) by 
P g (k) = b 2 P(k), motivated by the peak biasing paradigm (Kaiser 1984). The gravitational 
effects of dark matter lead to velocity fields which amplify the redshift space power spectrum 
in the linear regime. By measuring the distortions and imposing a stronger assumption 
5 g = b5, a parameter (3 = tt 6 /b can be constrained (Kaiser 1987). Recent attempts to 
measure these have resulted in a confusing picture, wherein the results depend strongly on 
the galaxy types (Hamilton 1997), scale (Dekel 1997) and surveys. A potential explanation 
for this state of confusion is that the linear galaxy biasing model is too simplistic. It has 
been proposed that galaxy formation may be non-local (Heyl et al. 1995). A general model 
assuming only that galaxy formation is a stochastic process has been proposed by Dekel and 
Lahav (1997). Unfortunately, the most general stochastic distribution introduces several 
new free functions, which make an unambiguous measurement of the matter distribution 
seemingly impossible. In principle, higher order statistics (Fry 1994) can be used to to 
break some of the degeneracies, but the application of the full theory is still in its infancy. 

Nevertheless, the fact that the correlation properties are distorted in redshift space 
tells us that gravitational interactions are at work. The challenge is to translate these 
measured distortions into conventional cosmological parameters. An added complication 
arises from the fact that the redshift distortions are more accurately measured in the 
quasi-linear regime where non-linear corrections are believed to be important (Hamilton 
1997). If we are given complete freedom to create galaxies with any distribution we please, 
and to allow these galaxies to have any form of correlation or anti-correlation with the dark 
matter density, can we still say something about the dark matter distribution? Can we 
invoke more observables to constrain the exhaustive class of galaxy distribution models? 
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In this paper we present the framework for a formalism which allows us to parameterize 
all galaxy distribution freedom into a simple unified picture. In section |2| we begin by 
constructing a hierarchy of approximations. The central limit theorem argues that the 
density field of both the galaxy distribution and of the matter distribution will tend to be 
Gaussian when smoothed on sufficiently large scales. Their correlation coeffient, however, 
remains a free parameter. For general random processes, a bivariate Gaussian distribution 
will have an arbitrary correlation coefficient — 1 < r < 1. The distribution has three free 
parameters: two variances and a cross-correlation coefficient. We show in section |3| how 
to determine all three using linear power spectrum distortions. The first deviation from 
Gaussianity is characterized by skewness. For a bivariate distribution, one has in general 
four skewness parameters. We show in section ^| how to reduce this number systematically. 
Once the power spectrum of the matter has been determined, the skewness of the matter 
can be computed using either second order perturbation theory (Hivon et al. 1995) or 
using N-body simulations. Skewness is also subject to redshift distortions, which can be 
measured. We show how these two additional observables allows one to constrain Q in 
section |3|. Individual galaxies may be a biased tracer of the local dark matter velocity field 
(Carlberg et al. 1990). But it is reasonable to expect that when smoothed on sufficiently 
large scales, the averaged velocity field will be unbiased. A geometric interpretation of the 
joint distribution function is given in section |4]. To assure accuracy, any parameters derived 
using this formalism must be checked against an N-body simulation. Such simulations also 
allow one to probe further into the non-linear regime. We describe in section |5| how to build 
mock galaxy catalogs which obey the formal expansion presented in this paper. 

2. Galaxy Distribution 

We will analyze the density field smoothed by a tophat window W R (r) on some scale 
R, which is zero for r > R and 3/AttR 3 elsewhere. The perturbation variable S = (p — p)/p 
is convolved to obtain the smoothed density field 5 R = J d 3 x.'5(x.')W(\x. — x'|). The galaxy 
density field will be described with a subscript g , and we similarly derive a smoothed galaxy 
density field 5 R . Only the galaxy field is directly observable. We do not have an exact 
theory about how and where the galaxies formed, nor how they relate to the distribution of 
the dark matter. The smoothed galaxy density field S R (x) could be a function of not only 
5 R (x.), but also of its derivative, and possibly long range influences. Even non-gravitational 
effects may play a role. Quasars exhibit a proximity effect, which may influence galaxy 
formation. If we consider all these variables hidden and unknown, we must prescribe the 
density of the galaxies S R as an effectively stochastic distribution relative to the dark matter 
5 R . 



-4- 



In the most general galaxy distribution, we need to specify the joint probability 
distribution function (PDF) P(5 R ,5 R ). At lowest order, we will approximate these to be 
Gaussians with some finite covariance 

P W , = ^exp(-^f! + M X-^). (D 

The parameters a, b, c are related to the traditional quantities by a change of variables: the 
variance of the matter a 2 = ({S R ) 2 ) = j J(5 R ) 2 Pd5 R d5 R = ajiac — b 2 ) and that of galaxies 
a 2 = c/(ac — b 2 ). The third free parameter is the covariance, given by the correlation 
coefficient r = (S R 5 R ) / a g a = b/^/ac. The traditional bias is defined as b\ = cr g /cr = Jc/a. 
The power spectrum of the galaxies is then related to the power spectrum of the dark 
matter by P g (k) = b\P(k). To simplify the notation, we will from now on implicitly assume 
that all fields are smoothed at scale R unless otherwise specified, and drop the superscript 
R on all variables. 

The analysis of this paper relies on Equation (|IJ) being a good approximation to the 
distribution of galaxies and dark matter when smoothed on large scales. A sufficient 
condition for this to hold is if the central limit theorem applies. Since we are considering 
the statistics of density fields averaged over some region, the central limit theorem will 
generally apply if the smoothing region is larger than the scale of non-locality of galaxy 
formation or non-linear gravitational clustering. There are some notable exceptions, 
however. 1. The dark matter power distribution may be non-Gaussian even when the 
fluctuations are linear. This might be true, for example, in topological defect theories 
of structure formation (Gooding et al. 1992) where the smoothed density field on any 
scale may be non-Gaussian. Such theories present great challenges to many attempts to 
measure cosmological parameters, including for example cosmic microwave background 
measurements (Pen, Seljak and Turok 1997). 2. The galaxy fluctuations could depend 
non-linearly on very large scale effects, for example external gravitational shear (van de 
Weygaert and Babul 1994). This could avoid the central limit theorem due to non- locality. 

Even if mild non-Gaussianity is present, we will argue below that the problem remains 
tractable. Conversely, there is every reason to believe that on small scales, the distributions 
of both the galaxy and the dark matter fields are significantly non-Gaussian. We will 
show in section |3| that mild non-Gaussianity actually allows us to break the degeneracy 
between the bias factor b\ and Q. We now turn to the next moment of the distribution, 
the skewness (5 3 )/a 3 . In principle one needs to specify four such independent moments 
(5 3 ~' l 5 l g ) for < i < 3. A coordinate transformation simplifies Equation (^ if we define 
y/2aS g = (u + v)/y/\ — r, \[2cb = (u — v)/y/\ — r, and w 2 = (1 — r)/(l + r). u is the 
variable with unit variance along the joint distribution of the galaxies and dark matter, 
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while v has variance w 2 and measures their mutual deviation. In this rotated frame, u and 
v are uncorrelated. In order to model all relevant terms, we apply a general Edgeworth 
expansion about the Gaussian (|I]). We recall (Kim and Strauss 1997) that the coefficients 
of the two first order terms u, v are zero since the mean is by definition 0. The three second 
moments are absorbed into the definitions of u, v. In principle one needs four third order 
Hermite polynomials to describe the joint distribution self-consistently at the next order. 
But we can reasonably assume that the galaxies are positively correlated with the dark 
matter distribution, so w < 1. In this case, the third order terms can be rank ordered in 
powers of w as u 3 , u 2 v, uv 2 , v 3 . As we will see below, it is necessary to retain the first 
two terms to model second order perturbation theory. We then neglect the last two terms 
because they depend on higher powers of the small parameter w, and disappear completely 
in the limit that biasing is deterministic r = 1. The Edgeworth expansion then gives us the 
truncated skew distribution 

D , s 1 + (u 3 - 3u)s + (u 2 - l)vb s /w 2 ( u 2 v 2 \ 

P(u,v) = exp . 2 

v ' ' Itxw f \ 2 2w 2 J w 

The two new coefficients are the joint skewness parameter s, and the second order bias b s 
which allows us to adjust the skewness of each distribution independently. The joint PDF 
will be discussed in more detail in section [| below. The Taylor expansion for general biasing 
introduces a quadratic bias at the same order as second order perturbation theory (Fry 
1994): 

5 g = f(5) = b x 8 + b 2 (5 2 - a 2 ) + 0(5 3 ). (3) 

We will show in section |] that in the stochastic notation, 62 = 1b s b\jo . 

We can now compute the basic relations needed for further calculations. All third order 
moments are uniquely defined 

(5 3 ) = a 3 (i±^) 3/2 (6 S -66 s ) 

/ 1 _l_ r \ 3 / 2 
(6 2 6 g ) = a 2 a g (^-) (6s-2b s ) 



2 

1 + r\ 3/2 
2 

1 + r\ 3/2 



^i) = (^+m s ) (4) 

For an initially Gaussian matter distribution with a power-law power spectrum 
P(k) = k n , the skewness factor of the evolved matter distribution is given by 
53 = (5 3 )/cr 4 = —(3 + n) + 34/7 from second order perturbation theory (Juskiewicz et al. 
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1993) which is exact for Q 
We obtain one relation 



1 and depends only very weakly on Q (Bouchet et al. 1992). 



1 + r 



3/2 



(6s — 6b s ) = a 



£-<3 + ») 



(5) 



2 

for the 5 unknowns a,a g ,r,s,b s . We will use five observational values to determine the 
remaining four galaxy distribution parameters, as well as the cosmological parameter Q. 
These determinations will rely on the comparison between redshift space and real space 
correlations. The real space correlation contains valuable information about the galaxy 
distribution itself, while the redshift space distortions are a consequence of the dynamics of 
the dark matter. 



3. Redshift space distortions 

We first consider the measurement of redshift space distortions of the variance or power 
spectrum in the presence of stochastic biasing. At first order, the two second order terms 
s, b s can be neglected. The redshift space density of galaxies is affected by their velocities 
through the Jacobian p g (x)dx = p g (x(z))(dx/ dz)dz and linear perturbation theory gives us 

S g (z) = 5 g (x) + 6(x)Q- 6 p 2 (6) 

where p = cos(6 l ) determines the angle between the wave vector k z and the line-of-sight 
(Kaiser 1987) and we have made the distant observer approximation. The power spectrum 
is the expectation value of the square of the Fourier transform of (|6|) and results in 

P g (k ZJ p) = P g (k)(l + /?V + 2r/V ) (7) 

where (3 = Q' 6 /b\. P g (k) is the undistorted power spectrum. The Legendre relation 

2/ + 1 r 1 



2 +1 Z" 1 

p i = ^~ J p^^Wp 



where Vi(p) is the /-th Legendre polynomial (Hamilton 1997) allows us to obtain moments 
of the angular dependence of (0). One can in principle measure both (3 and r by measuring 
the quadrupole distortion P 2 and the next order distortion P 4 separately. We can then solve 
for r and (3 using the following relations 

P 2 _ lr(3+l(3 2 



P 1 + |r/3 + 



P± Ml (9) 



P l + fr/3 + i/3 2 ' 
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With sufficiently large data sets, one can measure r and (3 as a function of wavelength fc. 
An alternative approach would be to measure P g from the angular correlation function, 
after which determination of the monopole Po/Pg and quadrupole Pz/Pg terms would be 
sufficient. Peacock (1997) compared P g derived from the APM angular power spectrum 
to Po to determine (3q = 0.4 ± 0.12 by setting r = 1. Allowing r to vary will increase the 
inferred value of (3 for all such measurements. In this case, the relation between the actual 
value of (3 for a given inferred value of (3$ is (3 = \J (3$ + 2(3 + r 2 — r. A similar increase in (3 
for a given inferred (3q using r = 1 holds for quadrupole measurements using equation @. 
If stochasticity has been neglected, all inferred values of (3 are only lower bounds. 

The skewness can be obtained from the bispectrum B g (ki, k 2 , k 3 ) = 
(5 s (ki)5 s (k 2 )5 g (k 3 )). Isotropy requires k x + k 2 + k 3 = 0, using which we can compute the 
third moment 

® = // d j2ny 2 B9{ku k2 ' ~ kl ~ ^W^W^W^k, + k 2 |). (10) 

One can then measure the net skewness by inverting the angular bispectrum to obtain the 
three dimensional bispectrum, in analogy to the power spectrum from APM (Baugh and 
Efstathiou 1994). The skewness of the smoothed galaxy field determines the last equation 
in (|j). The skewness of the dark matter field is determined by the variance allowing 
us to solve for both s and b s . Equations fl5[,|G|,|l0|) allowed us to solve for all parameters of 
the stochastic non-linear biased galaxy distribution model. The final goal is to break the 
degeneracy of (3 = Q s /b to determine Q and b independently. 

In redshift surveys, the measured skewness is already distorted by velocity distortions, 
but is nevertheless readily measurable (Kim and Strauss 1997). The second order 
perturbation theory calculation of the skewness has recently been completed for both 
redshift space and the real space (Juskiewicz et al. 1993, Hivon et al. 1995). Their 
basic result was that in the absence of biasing, the skewness factor S3 (see above) is 
weakly dependent on Q, and is very similar in real and redshift space in second order 
theory. In these calculations, Ss yZ and S3 typically differ by a few percent depending on 
Q . A similar sized change occurs when the bias parameter is changed. Fortunately, the 
observations already allow very accurate determinations, for example Kim and Strauss 
found 5*3 = 2.93 ± 0.09, where the errors are comparable to the expected effect. The actual 
redshift space skewness distortions quickly grow significantly larger than the second order 
predictions. For a self-similar power spectrum with n = — 1 and Gaussian filter o = 0.5 
Hivon et al. found S3 = 3.5 and Ss z = 2.9 using N-body simulations, which is a very 
significant effect and many times larger than current observational errors. Second order 
perturbation theory appears to systematically underestimate the redshift space skewness 
distortions. The real space skewness in simulations tends to be higher than perturbation 
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theory, while redshift space skewness tends to be lower in simulations. This trend suggests 
that direct N-body simulations are necessary to quantify the effect of redshift space 
distortions of skewness. We will examine this strategy in section [5] below. 

By measuring the skewness from solely angular correlation information (Gaztanaga and 
Bernadeau 1997) as well as from the redshift space distribution, we obtain two measures 
of skewness which can be compared to each other, from which one can solve for Q. One 
could also smooth the density field with an anisotropic window function, and compute the 
dependence of the skewness of the smoothed density field on the alignment of the window 
as was done for variance measurements by Bromley (1994). The formalism of Hivon et al. 
(1995) can then be modified using these anisotropic window functions. We will explore this 
approach further in section ^ below. Unfortunately, the second order perturbation theory 
redshift distorted skewness does not have a simple closed form expression, and a numerical 
triple integral must be evaluated for each specific set of choices of n, Q, W(k,f/,). Details 
of this procedure are given in Hivon et al. (1995). 



4. Coherent Limit 

A finite truncation of the Edgeworth expansion may result in a PDF which is not 
positive everywhere. For small corrections s < 1, 6 S < w 2 , the PDF in Equation @ 
remains positive in all regions where the amplitude is still large, and for practical purposes 
the negative probabilities do not have a significant effect. But it is possible to leave the 
regime of small corrections. When the corrections become large, the PDF becomes negative 
when it still has a significant amplitude. We can, however, absorb the coefficient of v into 
the exponent for small b s 

P(u, v) = 1 + ~ 3M) V" 2 / 2 X J? I (11) 

Equation (|TT| ) remains positive for all values of b s . The right term becomes a Dirac delta 
function in the limit w — ► 

P (u )V ) = l + (^-3^ e _ u2/2 x 6d[v _ ( u a _ i) 6s] . (12) 



We then reproduce Equation (§) to leading order for the choices b\ = ra/a g and 
b 2 = 2b s b 1 /a. 

It is instructive to understand the nature of the two skewness parameters s and b s 
graphically. In Figure [I] we show the joint PDF. The respective projections onto the galaxy 
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and dark matter PDF's is shown in Figure |2|. The projected PDF for the galaxies is 

i+gLJ^ f _£) (13) 

where s g = [(1 + r)/2] 3 / 2 (s + b s ) and that for the dark matter is the same with a change in 
sign of b s . For real surveys, one can measure the residuals of the galaxy PDF after fitting a 
third order Edgeworth expansion to determine the accuracy of the fit, with proper account 
for noise (Kim and Strauss 1997). The same can be done for the dark matter by utilizing 
N-body simulations described in the next section. 



5. Constructing Mock Catalogues 

Let us now examine how one can build galaxy catalogues from an N-body simulation 
consistent with this Edgeworth expansion. The purpose of this exercise will be to 
test specific models against catalogs in the non-linear regimes. We will show a sample 
construction of a galaxy density field smoothed by a window function which is consistent 
with the stochastic non-linear biasing described above. We can then compite the likelihood 
function for the cosmological parameters for any set of observations. The ultimate test 
would be to recover the same cosmological parameters and dark matter power spectrum 
using galaxy types which are known to be biased relative to each other. 

The first step is to calculate the bias function b(k) in Fourier space by comparing the 
angular correlation function of the survey to that of the simulation. Statistical homogeneity 
and isotropy require that the bias is a function of the magnitude of the wave-number only. 
We will take the density field of the simulation in Fourier space and produce the galaxy 
field by scaling to the bias 

^(k) = 6 1 (|k|)5 DM (k). (14) 

The mock density field is then convolved with the survey geometry and projected onto an 
angular power spectrum w(k). We solve for the bias function bi(k) by requiring the mock 
galaxy angular power spectrum to agree with the observed angular power spectrum. This 
procedure has used no velocity information. By repeating the simulation many times with 
different random seeds, we can obtain the full distribution of bi(k). 

We have three remaining parameters r, s, b s which must be solved for using velocity 
and skewness information. Since the skewness of the dark matter in the simulations is 
known, we have one constraint from Equation (f|), reducing them to two remaining degrees 
of freedom. We will use four observational quantities to constrain them: two moments of 
the redshift space variance distortion, and the skewness as well as its distortion. While 
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second order perturbation theory in principle allows us to solve for r, s, b s and fi, its 
validity quickly breaks down as one enters the non-linearly regime. We must use N-body 
simulations at this point to make quantitative comparison. The problem is now doubly 
overconstrained, allowing us to solve for two free simulation parameters, for example Q and 
the power spectrum shape parameter T by performing a sufficiently large number of N-body 
simulations (Hatton and Cole 1997). Velocity space distortions can be measured using 
anisotropic smoothing windows (Bromley 1994) W R (p). The window function proposed by 
Bromley symmetrizes the distribution and thus destroys skewness information. Consider 
instead an elliptical top-hat with a major-minor axis ratio of 2:1. 

We pick a characteristic smoothing scale R and smooth the observation on that scale. 
The trade-off occurs between picking large R, which smoothes over large volumes and 
results in distributions which are closer to Gaussian, or small R which results in a smaller 
cosmic variance and a stronger non-linear signal, but for which the first three orders of our 
the Edgeworth expansion may be a poor approximation for the true dark matter-galaxy 
joint distribution function. We first decohere the galaxy density field p g by adding an 
independent random Gaussian galaxy field p 1 ^ with identical power spectrum weighted by 
the correlation coefficient r: 

p' g = rp g + (1 - r)pf . (15) 

We have averaged the result of an N-body simulation with a random field with identical 
(non-linear) power spectrum. This maintains the shape of the power spectrum, but weakens 
the degree of correlation between galaxy and dark matter fields. It is no longer true that 
(S g \S) = bi5, but instead (5 g \5) = b\r5. Second order bias is added by feeding the field 
through a quadratic function 

ft = ^ + 7T H - prf - h2 A ■ ^ 

Next we distort into velocity space as follows: Each N-body particle mass is multiplied 
by p'g/pDM and projected with its velocity into a redshift coordinate system. The window 
function is applied in redshift space 



n i 

p R (z) = J p z (z')W R (\z - z'U)dV. (17) 

C(z) is the particle shape, which for Cloud-in-Cell mappings (Hockney and Eastwood 1980) 
is the same shape as the grid cell, n is the ratio of number of particles to the number 
of gridcells and is the scaled particle mass. Xi is the particle position, and vf is the 
line-of-sight component of the particle velocity which affects the radial redshift position. 
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We now compare the statistics with the observed sample. One computes the variance 
<r 2 (/i) = J(p R — p R ) 2 d 3 z and decomposes it into multipoles a 2 (p) ~ cr + o^Piip) + o" 4 P 4 (/i) 
as in Equation (^) and does the same with the skewness s^(p) = J(p R — p R ) 3 d 3 z where now 
S3 ~ s + S2P2(p)- These 02, 04, so, S2 are then compared with the values obtained from 
the surveys. A Monte-Carlo array of simulations provides the full likelihood distribution of 
these variables, allowing us to test consistency of each model with observations. 

This model of skew biasing allows us to discuss the systematic errors in the measurement 
of pairwise velocity disperions (Guzzo et al. 1997). Since pairwise galaxy velocities are 
measured in the non-linear regime, the inferred mean galaxy velocities can not be directly 
translated into mean dark matter pairwise velocity. Decoherence, and non-linear bias both 
introduce complex dependences in the conversion from galaxy velocity to dark matter 
velocity. Guzzo et al. (1997) showed that the one dimensional pairwise velocity varies by 
galaxy type from 345 to 865 km/sec. We must keep in mind that each galaxy type surely has 
different biasing and coherence properties. The pairwise velocities are typically measured at 
a separation of 1/h Mpc, where the density field is strongly non-linear, and the Edgeworth 
expansion may be a poor approximation to the actual joint galaxy-dark matter distribution. 
The mock catalog from N-body simulations described above effectively provides a handle to 
probe the dynamical properties of galaxies at larger separation, allowing us to separate the 
distribution properties of galaxies from the dynamical aspects of the dark matter. 

6. Conclusions 

The general stochastic galaxy biasing problem contains more free parameters than can 
easily be measured in any galaxy redshift survey. We have shown that using only linear 
perturbation theory we can determine two parameters, the correlation coefficient r and bias 
parameter (3 using the quadrupole and octupole distortions. This allows a reconstruction of 
the power spectrum P(k)Q' 6 as well as determination of two galaxy formation parameters. 
In the plausible scenario that galaxies correlate strongly with the matter distribution, only 
two free additional parameters s, b s need to be introduced to quantify the skewness. Second 
order perturbation theory provides one linear constraint, and observations of the skewness 
of galaxies determines the second. By measuring the redshift distortions of skewness we 
can in principle determine both the true underlying dark matter power spectrum and 
the density parameter Q independently. This picture has incorporated both stochastic 
correlation and second order non-linear bias. We have shown how to extend the Ansatz of 
the Edgeworth expansion to general problems without relying on linear or second order 
theory. In an N-body simulation the same approach can be applied to directly compare 
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specific models to observations. This also allows us to probe deeper into non-linear scales. 
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Fig. 1. — Joint Probability Distribution Function Contours. The parameters for this plot 
are a — 1, bias 6 = 2, correlation coefficient r = 0.8, skewness s = 0.05 and non-linear bias 
b s = 0.1. The solid line is the contour at half central probability, while the dotted lines are 
at 1/4 and 3/4. The axes are in units of the dark matter standard deviation. 
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Fig. 2. — The projection of Figure [I] onto the galaxies (solid) and dark matter (dashed), s 
parametrizes their common skewness, while each distribution's skewness is proportional to 
s ± b s . The units are in standard deviations of the dark matter. For our choice 6 = 2 the 
galaxies have twice the standard deviation of the dark matter. 



